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1. Scientific motivation for studying tropospheric OH 


3. Preliminary results for expected changes in OH TROP 


Sensitivity of OH TROP to 
climate change 

• IPCC (2001) predicted future OH would fall due to rising CH 4 

• Best knowledge of global tropospheric OH (OH TROP ) comes 
from analysis of CH 3 CCI 3 observations 

• The CH 4 lifetime inferred from OH TROP based on CH 3 CCI 3 is 8.9 
years; yet many sources use a CH 4 lifetime of ~12 years, which 
was found by modeling studies described in IPCC (2001) 

• The 12 year CH 4 lifetime is central to subsequent IPCC reports 
and is called the "perturbation lifetime" 

• Our work is motivated by understanding: 

- what factors other than rising CH 4 will affect OH TROP 

- veracity of notion that future OH TROP will decline 

Importance of OH TROP 

• OH serves as the main sink for many species: CO, CH 4 , 

CH 3 CCI 3 , all HCFCs and HFCs, biogenic VOCs (isoprene and 
terpene), anthropogenic VOCs (formaldehyde, benzene, 
toluene, etc.) 

• HO x is integral to 0 3 formation and destruction pathways; 
alone, HO x depletes 0 3 , but with NO x and VOCs, HO x creates 0 3 


Methane Oxidation and OH 

• Under low NO x , the oxidation of CH 4 destroys HO x 

• Under high NO x , the oxidation of CH 4 produces HO x 

• Loss of methylperoxy radical CH 3 0 2 with NO or H0 2 is critical 

• Future levels of CH 4 and NO x are highly uncertain 

- CH 4 between now and 2100 varies dramatically among 
RCP scenarios 

- NO x between now and 2100 will depend on whether the 
developing world implements selective catalytic reduction 
on coal power plants and catalytic converters on cars 

• Future OH TROP will also depend on overhead column 0 3 , 
local humidity, and biogenic emission of VOCs, etc. 

• Given this complexity, we expect CCMs to project 

a wide range of values for OH TROP during the rest 
of this century 2000 


NO x emissions from industry, 
power generation, traffic, 
domestic heating, and biomass 
burning used as input for prior 
CCM calculations, for years 2000 
(38.0 TgN/yr total) and 2030 
(67.6 Tg N/yr total). From 
Eyring et al., ACP, 2007. 




Previous Work on OH TROP Trends 


• Montzka etal., 2011 found that OH TROP does not vary interannually (from 1997 to present) 

• Prior studies by Prinn et al., 2001 and Bousquet et al., 2005 suggest large interannual variability in 
OH TROP (1985-2000) 

• We suggest the OH TROP behavior in all three studies may be physically possible, based on our 
preliminary results 


Current Conclusions from OH TROP Analysis 

1 Effect of CH 4 on OH TROP is taken from IPCC (2001), Section 4.2. 1.1, which states that "the feed- 
back of CH 4 on tropospheric OH" found using contemporary chemical transport models is -0.32% 
for every 1% increase in CH 4 (red line, middle panel of figure to right) 

1 Primary effect of overhead 0 3 is rise in OH TROP following the 1991 eruption of Mount 
Pinatubo (green dashed line) due to enhanced removal of stratospheric 0 3 by volcanic aerosol 

1 Rising H 2 0 from MERRA and AIRS increases OH TROP (orange lines) by an amount comparable 
to the decrease expected from rising CH 4 

Widening of the tropics contributes an increase in OH TROP (purple line) roughly equal to that due 
to rising H 2 0, leading us to view this as a significant factor controlling OH TROP 

1 Overall expected change in OH TROP (bottom panel) shows higher level of interannual variability 
prior to ~1999 and lower variability thereafter 
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2. Methods for estimating changes in OH TROP 


4. Future plans for refining estimates of AOH TROP 


Model OH 

Initial modern-day conditions for OH 
are taken from a time-slice run of 
GEOS CCM using 2005 emissions 

Monthly mean mixing ratios of OH 
and related species are provided on 
a 144 longitude, 91 latitude, 72 
pressure level grid 
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Changes in OH due to H 2 0, 
overhead 0 3 and tropical widening 
are found from analysis of 
initial OH field 



MERRA & AIRS H 2 0 

• Specific humidity files from the NASA 
Modern-Era Retrospective analysis for 
Research and Applications (MERRA) 
starting prior to 1985, & 

• H 2 0 mixing ratio files from the NASA 
Atmospheric Infrared Sounder (AIRS) 
starting in 2002 were used to establish 
trends in H 2 0 by latitude 

• OH TROP is assumed to follow the square 
root of the change in H 2 0 using a 
steady-state assumption 

• We plan to refine the effect of changing 
H 2 0 on OH TROP by examining reaction 
rates from archived runs of GEOS CCM 


Overhead 0 3 

• Total column 0 3 trends were obtained 
from the NASA merged 0 3 data set, 
consisting of measurements from SBUV, 
TOMS, and Aura OMI instruments 

• We then use our photolysis code to 
estimate the impact on J(0 3 ) O^D) of 
decreasing initial GEOS CCM overhead 
0 3 columns by amount suggested by 
the NASA product 

• OH TROP is assumed to change by the 

square root of J(0 3 ) 0( 3 D) 



Tropical Widening 

Widespread evidence for expansion 
of the latitudinal circulation system 
(e.g. Seidel etal., 2008) 
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GCMs forced by GHGs cannot 
reproduce extent of tropical 
widening (Johanson and Fu, 2009) 

Allen et al. (2012) have suggested 
tropospheric 0 3 and BC are 
responsible 

We simulate tropical expansion by 
increasing OH near the tropical 
boundary by a factor representative 
of 2 “widening/decade (globally) 


Improve Estimate of 
d(OH TROP ) / d(H 2 0) 

Reaction rates from recent runs of GEOS CCM 
are archived for reactions such as: 

HgO + OOD) 20H 

Using these reaction rates we will determine 
the proportion of OH that is produced via 
reaction with H 2 0 

The determined scaling factor would be used to 
calculate a new AOH TROP based on the H 2 0 
trends 

Estimate time- and pressure-varying values of 
d(OH TROP ) / d(H 2 0) 

Evaluate discrepancies between MERRA and 
AIRS H 2 0 trends 


• Calculate uncertainties in AIRS and MERRA 
H 2 0 and NASA 0 3 product 


Evaluate CH 4 /OH Feedback 

■ We will use a box model (details below) to probe 
relationship between CH 4 and OH TROP and its 
dependence on NO x 

Box Model 

We will use the Dynamically Simple Model of 
Atmospheric Chemical Complexity (DSMACC) box 
model (provided by Mathew Evans) with chemical 
mechanistic information taken from the Master 
Chemical Mechanism, MCM v3.2 (Jenkin et al., 

1997; Saunders et al., 2003) via website http:// 
mcm.leeds.ac.uk/MCM. Using the MCM, we have 
the capability to consider: 

• the degradation of up to 135 VOCs 

• the influence of J(0 1 D) and J(N0 2 ) on OH 
chemistry 

• the impact of NO x on methane oxidation 


Evaluate standard deviation in average fraction 
of OH production occurring via H 2 0 + 0( 1 D) 

Estimate uncertainty in the box model 
evaluation of d(OH TR ° p ) / d(CH 4 ) 


Assessing OH TROP in CCMs 

• Through our involvement with the IGAC / SPARC 
Chemistry-Climate Model Initiative, we have 
requested: 

1. hourly, instantaneous output from participating 
CCMs 1 day/season, 1 year/decade 

2. archival of all species, reaction rates, J-values, 
and physical parameters relevant to OH 
chemistry 

3. this output for both the REF-C1 (hindcast) and 
REF-C2 (future) runs 

• We plan to assess the causes of differences 
between OH in the CCMs 

• Use of the box model enables us to distinguish 
between OH differences due to chemical 
mechanism and those due to differences in OH 
precursors 

• We can also use this output to predict future trends 
in OH TROP , based on CH 4 , H 2 0, and overhead 0 3 from 
the future CCM runs 


Propagate Uncertainties 







